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Abstract 

We present an ab-initio formalism for the calculation of transport properties in compositionally 
disordered systems within the framework of the Korringa-Kohn-Rostoker non-local coherent po- 
tential approximation. Our formalism is based upon the single-particle Kubo-Greenwood linear 
response and provides a natural means of incorporating the effects of short-range order upon the 
transport properties. We demonstrate the efficacy of the formalism by examining the effects of 
short-range order and clustering upon the transport properties of disordered AgPd and CuZn 
alloys. 

PACS numbers: 72.10-d, 71.20.Be, 72.15.Eb 
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I. INTRODUCTION 



The Korringa-Kohn-Rostoker coherent potential approximation (KKR-CPA) 1 represents 
an extremely successful one-electron theory capable of describing the properties of many 
compositionally disordered alloy systems. In particular, in combination with density func- 
tional theory (DFT)^ 1 ^ and the Korringa-Kohn-Rostoker method of band theory***, it 
provides a fully first-principles description of such systems. A long history of successful 
applications^^^^ 1 ^ attests to the utility and accuracy to which the method is capable. 

Of particular relevance to our discussion here, the KKR-CPA has been used to calculate 
the transport properties of such alloys. Historically, the early work concerning this topic 
involved the use of a Boltzmann equation^^^, and although these works demonstrated 
remarkable agreement with experiment, they did suffer from two notable defects: namely, 
the requirement that well-defined energy bands exist within the alloy, and the neglect of 
vertex, or "scattering-in" , terms. The former defect limits the application of such a theory 
to weak-scattering alloys only, while the latter could be expected to lead to significant error 
in systems such as those where appreciable s-p or s-d scattering manifests itself. 

Velicky- 8 developed a CPA theory for transport using the Kubo-GreenwoocU*^ formalism 
as his starting point, which was applied to a two-level tight-binding Hamiltonian. Although 
capable of yielding vertex corrections, this approach suffered from difficulties when applied 
to realistic systems. For example, the necessity of assuming that the wavefunctions are 
identical on all lattice sites, irrespective of the occupying atomic species. The seminal work 
of Butler—, resolved these issues, as he developed a KKR-CPA theory based upon the Kubo- 
Greenwood linear response formalism^^. As a multiple-scattering based approach, this did 
not require the existence of well-defined energy bands, and further, Butler demonstrated 
how the vertex corrections arose quite naturally within his formalism. Although the formal 
developments of this work followed Velicky's quite closely, the use of a realistic single-electron 
muffin-tin Hamiltonian allowed connection with first-principles methods to be made. The 
method has been applied with success to a range of alloy systems^ 1 ^ and it has also been 
successfully extended to the relativistic regime^. 

Of course, all of these calculations suffer from the main drawback of the CPA; namely 
that as a single-site mean-field theory, it is incapable of incorporating the effects of local- 
environment fluctuations in the alloy crystal potential. The CPA therefore explicitly ignores 
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the effects of such short-range order (SRO) effects upon the physics of disordered alloys. 
As discussed by Gonis 25 , these statistical fluctuations can be important. Whilst in general 
the presence of SRO is likely to diminish the resistivity there are examples, the so-called 
"Komplex" K-state alloys 26 , where the onset of SRO is accompanied by an increase in the 
alloy residual resistivity. 

Recently there have been some successful attempts at calculating the effects of SRO upon 
the electronic structure of disordered alloys and they provide the means to study transport 
properties. Mookerjee and Prasad 2 ^* 2 ^ 2 ^, using a TB-LMTO method 3 ^ in conjunction with 
an augmented space formalism^ 3 - 2 - and real space recursion method 33 described SRO effects 
on alloy electronic densities of states and related quantities whilst Saha et a/.— have obtained 
spectral functions within the same framework. Recently Tarafder et al.— have developed 
a formalism for the optical conductivity and reflectivity from the same basis and used it 
to study copper-zinc alloys. To date however it has not been possible to incorporate this 
technique fully within electronic density functional theory. The recent work of Rowlands et 



a p£i2L2§ 7 along with Biava et a/.— ^ concern a development which is not restricted in this 



way. References 
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functional theory to provide a first principles description of disordered alloys as demon- 



strated in reference 
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This nonlocal CPA (NLCPA) theory is based upon reciprocal 
space coarse-graining ideas introduced by Jarrell and Krishnamurthy- 2 originating from the 
dynamical cluster approximation (DCA)^ 3 - 1 ^ 1 ^. The KKR-NLCPA— ^ introduces an effec- 
tive (translationally invariant) disorder term SG which represents an effective propagator 
that accounts for all nonlocal scattering correlations on the electronic propagation due to 
disorder configurations and modifies the structure constants accordingly. By coarse-graining 
reciprocal space, one naturally introduces real space periodically repeating clusters. As such, 
the NLCPA maps an effective lattice problem to that of an impurity cluster embedded in a 
self-consistently determined effective medium, and thus yields a cluster generalisation of the 
KKR-CPA that includes non-local correlations up to the range of the cluster size. Unlike 
other cluster approaches, such as the molecular CPA (MCPA) 46 it is fully translationally 
invariant, that is, the effective medium has the site-to-site translational invariance of the 
underlying lattice. It is also computationally tractable, largely on account of the reciprocal 
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space coarse-graining procedure employed. 

Thus far, the NLCPA has been employed to investigate the effects of SRO upon the elec- 
tronic structure of a range of realistic systems, using both the muffin-tin Hamiltonian2I^, 
and the tight-binding approach^. Given its ability to address such issues successfully and 
its proven incorporation into an electronic density functional theory^ (DFT) for disordered 
systems with SRO it makes sense then to extend it to the calculation of transport proper- 
ties. To this end, by invoking time dependent DFT — within the adiabatic approximation, 
in this paper we present a formalism for the determination of the residual resistivity in the 
KKR-NLCPA, and explicitly demonstrate the efficacy of the method through application 
to several realistic alloy systems. Our theoretical formalism is a careful generalisation of 
that of Butler— and where appropriate we omit the steps in the derivation which can be 
straightforwardly obtained from this paper. Our paper is structured as follows: The next 
section gives a short overview of the conductivity tensor and then the transport coefficients 
available from the Kubo-Greenwood formalism. This is followed by a section containing the 
salient points of the KKR-NLCPA formalism including its use in implementing the density 
functional theory. We then develop our theory for the conductivity of disordered systems 
with short-range order which includes the treatment of 'vertex corrections'. The implemen- 
tation strategy is outlined before calculations for the effects of SRO on the resistivity of 
both b.c.c. CuZn and f.c.c. AgPd are presented. 

II. THE CONDUCTIVITY; KUBO-GREENWOOD LINEAR RESPONSE 

The Kubo-Greenwood^ 1 ^ linear response formalism states that, for disordered system, 
the symmetric part of the conductivity tensor has coefficients C which can be determined 
from the evaluation of an expression of the form 



where G is a single-particle Green's function, which is dependent upon the details of the 
effective one-electron potential. 0\ and O2 are operators, and the angled brackets denote a 
configuration average over the distribution of the potentials. 

To determine the d.c. conductivity, we consider the following expression 



C = Tr < 1 G0 2 G >, 



(1) 




mn 
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where J mn;A1 = {m\J fJ \n) denotes the matrix element of the current operator in the /zth 
spatial direction, which is given by 

6 

J» = -ih—-^—, (3) 
m or^ 

with \m) and |n) denoting eigenfunctions of a particular configuration of the disordered 
system. Here, N is the number of atoms, and Q is the volume per atom. Ep is the Fermi 
energy. 

Within the KKR approach, the electronic structure of the alloy is expressed in terms of 
the single-particle Green's function, rather than in terms of eigenstates and eigenvalues of 
the Hamiltonian; we can introduce the Green's function simply by using the identity 50 

-ttV \n)(n\5(E - E n ) = \imlmG(E + irj), (4) 

while the awkward imaginary part of the Green's function may be removed by writing 

-2iirJ2\n)(n\5(E - E n ) = ]xm[G(E + irj) — G(E - irj)]. (5) 

Inserting this into Eq. (T5]), yields the following for the conductivity 

<V = l/4lim[d^{E + ,E + ) -a^(E~,E + ) - 

d^ u (E + ,E-)+d, u (E-,E-)] (6) 

where we define the complex energies as 

E + = E F + iT], E~ = E F -ir), 7] -> (7) 

and 

~o»M, Z2) = |^Tr < J fl G(z 1 )J v G(z 2 ) >, (8) 

where z\ and z 2 are each either E + or E~ . 

For non-overlapping effective single electron potentials the Hamiltonian takes the form 

H=-^V*+Y,v 1 {T-ll i ), (9) 

where the atomic positions R» are fixed, and form a regular lattice. The potentials u 7 (r$) 
vary from site to site (r^ = r — Rj) and 7 is a configuration label. 

Within the multiple-scattering theory, the single-particle Green's function for a given 
configuration 7 can be written as^i 



G 7 (£,r 4 ,r;.) = 2m/hH £ Zi^E^r^Z^E^ - £ ^ 7 (r<) J^ 7 (£, r>)^ ,(10) 

V A, A' A J 

where t*\ a/ is the scattering path operator (SPO) describing propagation between sites i 
and j in configuration 7, Z l A ^(E, r») is the regular solution to the Schrodinger/Dirac equation 
in the cell surrounding the atom i, and J A7 (i?,r^) represents the irregular solution within 
the same cell (note that there should be no confusion with the current matrix elements here). 
A encapsulates the appropriate angular momentum quantum numbers.—" 2 ^ 

In calculating the conductivity, the second term in the Green's function expression (Eq.( 
[TO]) ) is real and may be omitted, when calculated for a real potential at a real energy^. 
Thus the conductivity may be written 



Am 



with 



E E <JZ 



i,j Ai,A 2 ,A3,A 4 



^AA>( Z ' z ) — 



ieh 



A 1 A 2 ,^K Z 2, Zl) T r,A 2 A a (zi)JA]A 4 ,u( Zl i Z 2)t^;A 4 A 1 (^2) > 



(11) 

(12) 



^ r «^A,7( r «5 Z )~Q^r^A',-y( r iy z ) 



771 Jcelk 

where celU defines the region surrounding the site i. 

We now need to consider how to carry out the averaging over configurations implicit in 
Eq. ljTTI) . Butler— showed in detail how to use the CPA to accomplish this. The single site 
nature of this effective medium theory means, however, that the potentials on the different 
lattices could only be treated as statistically independent. We will show how to carry out the 
averaging using the NLCPA whereby short-ranged correlations can be naturally included. 
To this end in the next section we summarise briefly the key aspects we need from the 
KKR-NLCPA together with its incorporation into electronic density functional theory. Full 



details can be found in references 
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III. THE KKR-NLCPA AND ELECTRONIC DENSITY FUNCTIONAL THEORY 

In general the scattering path operator (SPO) between two sites i and j for an electron, 
moving through an effective medium so that it mimicks the average motion in a disordered 
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system, is given by 



T 



tSij +E*( G ( R i - R fc) + 5G ( R i ~ R k))f k3 . (13) 

Here, all quantities are matrices in angular momentum space and the indices i,j run over 
all sites in the lattice. The G(Rj — Rfc)'s are structure constants. The effective medium is 
specified by single site t-matrices t and effective structure constant corrections 5G(TLi — Rj). 
The effective medium must be translationally invariant so that f u is given in terms of a 
Brillouin zone integral. 



f dklr 1 - G(k) - ^(k)]-^*-^-^^. (14) 

IIbz J^bz 

In order to establish a tractable procedure for determining the effective medium the NLCPA 
draws its chief idea from the dynamical cluster approximation (DCA) for interacting elec- 
tron systems^ 1 ^. This is a coarse graining of SG consistently in real and reciprocal space 
and a mapping to a self-consistently embedded impurity cluster problem with appropriate 
boundary conditions imposed^. The full translational symmetry of the underlying lattice 
is preserved. The size of the cluster sets the range of correlations that can be included. 
The lattice is divided into 'tiles' centred on a superlattice vectors Rc and each contains N c 
sites at positions Rc + R/, / = 1, • • • N c . The Brillouin zone is also broken into N c tiles, of 
volume Q t — &>bz/N c , centred on the cluster momenta K n , n = 1, • ■ ■ N c and the R/s and 
K n satisfy the following equation: 

1 £ eJ K„.(R,-R j)=5/j _ (15) 

5G(k) is coarse-grained so that it has the average value 5G(K n ) in a tile centred on K n and in 
real space, <fG(Rj - Rj) = (l/N c ) £ Kn 5G(K„)e iK " ( R ^ R ^ with 5G(K n ) = Ej^G(R 7 - 

p_ ^ e -iK„.(R x -Rj)_ 

The SPO is coarse-grained, 

Nc 

t-<BZ 

appropriate to the reciprocal space tile of volume Q t and in real space for multiple scattering 
starting and ending on cluster sites I and J respectively. 



f(K n ) = — 1- / dk.lt- 1 - G(k + K n ) - SGiKn)]- 1 (16) 



f IJ = J- 2( / ^[r 1 - G(k + K n ) - dGiK^}- 1 )^^-^. (17) 



Note how e* k '( R/-Rj ) is taken to be ~ 1 as the coarse graining is applied 43 . The final step 
is to find the SPO for an impurity cluster, describing a particular configuration 7c of 
atoms which is embedded into the NLCPA medium. By demanding that the average is equal 
to the SPO of the NLCPA enables the effective t-matrix and structure constant corrections 
to be determined, i.e., 

E*V£ = *"- (18) 



1C 



Short-range order, SRO, can be inc 
demonstrated in, for example, refs. 



uded by choosing the probabilities P 7C appropriately as 
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In ref. 
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in a generalisation of the work of Johnson et al. 4 ' 5 , it is described how to 
specify a configurationally averaged electronic Grand potential, Cl in terms of KKR-NLCPA 
quantities and charge densities p 7c (rj) and one-electron potentials, v lc (ri) different for each 
cluster configuration. The functional minimisation of with respect to the charge densities, 
p 7c (r/)'s, determines the total energy of the system and requires the p 7c (r/)'s and v lc {vi) 's 
to be found self-consistently. Rowlands et al— applied this DFT to investigate how the total 
energy, charge densities and densities of states are affected by SRO. Tulip et al.— showed 
how further information can be found about the effects of SRO on the electronic structure 
by formulating and calculating the Bloch spectral function at the cluster momenta and 
averaged over tiles whilst Batt and Rowlands 4 ^ explained how the spectral function at any 
point in the Brillouin zone can be found. In the following we build on these developments 
and describe the theory for a two-particle correlation function of disordered system with 
SRO. The particular example is to the d.c. conductivity. 



IV. ANALYTIC CONFIGURATION AVERAGING OF THE CONDUCTIVITY 
USING THE KKR-NLCPA 

From the above it is clear that the KKR-NLCPA should enable an analytical configu- 
rational average of the conductivity to be carried out. To do so, some care needs to be 
exercised. There are two distinct cases that we must consider in Eq.flTTJ: (i) where the two 
sites under consideration, % and j, lie within the same NLCPA cluster, and (ii) when they 
lie in two different clusters, in which case the occupancies of the sites will be statistically 
independent, as the two distinct clusters will be statistically independent. This is a natural 
generalisation of Butler's work—, where he distinguishes between the two cases of i — j and 
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Hereon we use lower case letters to denote general sites in the lattice, ■ ■ ■; upper case 
C denotes tiles containing the clusters; and upper case letters, with the exception of C, 
denote sites within clusters. So for a site at position R; we use Rj = Rc + R/ and I, J, - ■ ■ 
label sites within tile C, I', </',••• sites within tile C etc. 

We accordingly write 

&ia,(zi, z 2 ) = a (z u z 2 ) + a 1 (z u z 2 ) (19) 



Z 2) = --SiF, 12 12 < J AiA 2 ;^2, Z^rl^iz^jl^izu Z 2 )t^ Ai (z 2 ) > 

irn \i j e c Ai,a 2 ,a 3 ,A4 

(20) 

and 

^(^1^2) = 12 12 12 < 4 1 A 2 ^(^2,^l)rA;X +J '(^l) J AU4^( Z l' Z 2)rf 4 X' / '' / (^2) > 

Tin \L c'^G J'eC h x MMM 

(21) 

where <7° includes sites J within the same NLCPA cluster (denoted C) as our reference site 
J, and a 1 includes all sites lying outside this cluster. Note that in writing these equations, 
we have utilised the translational invariance of the averaged system to remove the second 
sum appearing in Eq. flTT]) . 

We now introduce response functions such that we can write 

3U*i, ^) = -^12 P-yo 12 JaX-A^ *i)k'£S£>>&> z *) ( 22 ) 
Klg$M> *>) = 12 12 < Tli A M)J J A:z; u (zu ^^(^ >c, 7c . (23) 

Jec A 3 , A 4 

K^ 2 \° u (zi, z 2 ) involves an average over all configurations except the configuration is fixed 
in cluster C to be 7c- The single site quantity J^^- J,^ z i) is se t U P depending on what 
kind of element occupies site I and the one-electron potential i> 7c (rj) that is dependent on 
the configuration 7c 1 . 

Similarly, for the inter-cluster contributions to the conductivity, we can introduce the 
following 
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^(21,22) — £ ^2 Pr rc^2 P 7o' £ ^a'X;/^ 22 ' z l) L A2,Ai;^ ' 7C "( z 1> z l) ( 24 ) 

Tin \L c ,^ c lc 7(JI Ai,A 2 

and 

T I,C,jc;C',-y c , i \ ^ I,C'+J' / \ rJ'ac 1 \ C+J',1, \ . / OK \ 

L A 2 ,A i;i , (21, 2 2 ) = £ Z^ <T A 2 A 3 (^lJ^AsA^^b^jT^ (Zz) >C,7c;C',7 C ' ( 25 ) 

J'eC" A 3 ,A 4 

where the notation is similar to before but the average now fixes cluster C to be occupied 
by configuration jc an d cluster C to be loaded with configuration 7^. 

In order to evaluate the ensemble averages contained in Eqs. (f2"2l) - (j25i) it is helpful to 
express the SPO, r, for a particular configuration as the SPO in the NLCPA medium plus 
corrections. Again this is the direct generalisation of Butler's approach^! We can write 

T H = f ij + ^2 f lk T kl f lj (26) 

k,l 

where the effective medium path SPO is denoted f as before and T is the total scattering 
matrix relevant to the specific configuration. The double summation is taken over all lattice 
sites, t 4 - 7 satisfies the following equation 

Y(tk%k - G(Ri - R k ))r kj = (27) 

k 

where all quantities are matrices in angular momentum space. We consider fluctuations 
about the NLCPA medium to obtain 

£ (W - r 1 )^ + 6G(R t - R fc ) + r 1 - 5G(Ri - R fc ) - G(R t - R*))t*' = 5 tJ (28) 

which can be re-arranged to yield 

T H = f ij _ T ik Am kl T lj . (29) 

k,l 

with 

Am kl = (t^ 1 - r l )5 kl - 5G(R k - R) (30) 

and have used the fact that the effective medium SPO may be written as the inverse of the 
matrix with elements (t _1 5ij — 5G(R.i — Rj) — G(R» — Rj)). 
We can thus write down 

J2T kl f lj = ^Am¥. (31) 
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If we now substitute for r' J using Eq.( [26]), label the sites according to clusters, C and sites 
within those clusters (uppercase letters), we obtain 

rpC+K,C'+L' _ x KLfi Gci _|_ y^ x KM~C+M,C"+N"jiC"+N",C'+L' ^ 2 ) 

C"^C MN" 

where we have introduced the matrix x associated with a single cluster of sites, given by 

x IJ = fl + Amf) Aw*- 7 . (33) 



A" 



JA 



Our results here are a direct cluster generalisation of Butler's^. Note also that these 
results are consistent with Hwang et al's cluster CPA conductivity formalism^ (although 
that is phrased in terms of t-matrices, rather than the x-matrix that we use in this work). 
Further, the special case of a single-site cluster recovers the more familiar CPA results. 

The NLCPA amounts to writing 



, rpC+K,C'+L' , , KL . r i \^ \^ ^ „AAf . ~C+M,C"+N" , T <7"+iV",C"+L' . 

< -t >vlcpa— < a; > oc,c+ 2^ < x >r <J >nlcpa 

C"^C MN" 

(34) 

and if we choose that < x IJ >= 0, which is another way of expressing the NLCPA ansatz 
(Eq. (IT5|) ) then < T >atlcpa= 0, and we obtain < r >nlcpa= t, (Eq. (ilSj) ). Of course, in 
writing this, we have made the approximation that 

< x KM~C+M,C"+N"rpC"+N",C'+L' > ~ < ^AA/ > ~C+M,C"+N" < rpC"+N" ,C'+L' > 

which is analogous to the usual CPA-type averaging approximation. 

Using Eq. (|26|) and closely following a cluster generalisation of Butler's analogous deriva- 
tion, which refers to fluctuations about the single site CPA medium, we find (suppressing 
angular momentum labelling) 

M,N 

with 

^ = E^*^ J + E E ^ c " +x "r^ LV f c "^"^ (37) 

AA C"^CK",L" 

and 

M,A" 



L v U IM Lj MN;v U NI \ 6 °) 
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where 

;C*,7 C ;C",7 C , _ \ ^ ~M,C'+K' rlc' ^.C'+L 
^MN-v — ' J K'V-y 



K'V 



+ E E ( fM,c +K >N )- ( 39 ) 

C"+(C,C>) K",L" 

In Eqs. (J3SD and (J3BD we use the NLCPA projector D(D ] ) which is D = (1 + Amf)" 1 
((1 +fAm)) found in Eqs. (1301) and (I33p . The NLCPA ansatz can be re- written in terms of 
them, i.e. X) 7C P-y C D lc = 1- We have also defined the current quantities 

J KL;u — Z^ LJ KN J u U NL- l 4U J 



N 



in Eqs.f l3"T|) and fl3"§|) . Finally we have introduced vertex functions r^',7^» L „. u , ^c'^'k" l'^-I 
which are the NLCPA analogues of the vertex functions derived by Butler 21 . We now show 
how to calculate these quantities. 



V. VERTEX FUNCTIONS AND THE CONDUCTIVITY 

Our two vertex functions are slightly different. ^c'" C k" v-v which appears in Eq.([3j 
for the intra-cluster component of the conductivity, concerns the connection between the 
configurational occupation in one cluster, C, with that of another, C", whereas the vertex 
function ^c'^jc'^u f° r the inter-cluster contribution relates the contents of cluster C" with 
that in two others, the reference one, C and another C . To facilitate the derivation of a 



21 



and 



closed set of equations, we introduce an approximation, analogous to that in ref. 
assume that the dependence of the latter on the contents of cluster C may be neglected. 
Thus 

r C,7c;C",7 c / _ F C",7 C , 

1 C",K",L";v — 1 C",K",L";V V±"U 

leading to L%$£' ,7c ' = in Eq.flM])- Using Eq. (J2ED, the NLCPA condition, Eq.flTg} 

(or its equivalent renditions) and tracking the steps in ref.— we obtain 



F C",7c . V <T T IK >~fC~K,C'+M> p c > ~C'+N',L LJ,-y c 

1 C,IJ;u — <~ x T J M'N';u T x ^C'^ c 

K,L,M',N> 

+ E E <a^f**' + ^ 

C"^(C,C) K,L,M",N" 



12 



which, if we compare with Eq. (l39l) . allows us to write the vertex function in terms of the 
response function 

r?SZ = E < x IK ™L C k %'x lj ™ > c , acl • (43) 

K,L 

This yields a closed set of equations for the conductivity. We may now write Eqs. ( |22l) 
and (JMD as 

A 7 Y) 

**) = —f??w E p 7c E E JJimm-A^ *i)K?jiX*>M> **) ( 44 ) 



and 



^IM,^) = - j! E E E E p 7c p 7c x 

7rrt iijv c i e c C'^C J'eC'icya' 

EE 

(zi,Z2). (45) 

Ai,A 2 /,J 

The response functions that determine the conductivity are given by: 

Klj!kiM\v = EE ^A^3^i,A3A4;v^A 4 J A2 + (46) 

A",L A 3 ,A 4 

E\r- ^ >.I,C'+K' ~i c ,,K'M' jC,~f C lc>> N ' L ' ^ *C'+L',J 

/ , Z^A 3 ,A 4i A 5 ,A 6 r AiA 3 < X A 3 A 4 ^M'N'MAs^AsAe > C,-yc r A 6 A 2 
CVC K',M',N',L' 



^IJMMW ~ AiA 3 J A-'L',A 3 ,A4^ r A 4 A2 + I 4 ' J 

A'',L' A 3 ,A 4 

EV^ V- ~/,C*"+A"' i c ,,,K"M" fC'tfc, 7c „,V",L" ~C"+L",J 

Z^A 3 ,A4 ! A 5 ,A 6 r A!A 3 < X A 3 A 4 iv M",Ar",A 4 ,A 5 ;i/ a; A5A 6 > C',lc> T A 6 A 2 

C"^(C,C) K",L",M",N" 



VI. SOLUTION OF TRANSPORT EQUATION 



The conductivity is determined by solution of the transport equation, Eq. fT47|) . We first 
see that in the key Eqs. (1511) and (H5|) we require the response function L^j?^' averaged 
over configurations jc 1 t na t can be assigned to a cluster C in a NLCPA tile located at 
position Re-/ and also that we require this to be summed over all clusters C . We rewrite 
Eq. dUD, omitting the quantum numbers, Ai,A 2 etc. for brevity, and denote the averages 
E 7a , Py cl Lfj'^ c ' , E 7c , ; P-i c ,Jk C '',l';u as an d Jk>,l>;u respectively as well as summing over 
clusters C 

L IJ;u = 1^ 1^ T Jk',L';uT ~ ^ T ' J K ',L';uT + 

C C K\L' K',L' 
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EE E ^^"vT^^^L^^ (48) 

C C"^C K",L",M",N" 

where w K ' ,L > M '< N =< x K " ,L " x M " ,N " > and Lfj. u is defined as zero. 

We now write the SPOs in terms of their lattice Fourier transforms, i.e., 

jJ,c>+J> = _L £ f dk^K n] k)e l ^ + V< ni ~ R ci-Kj>) (49) 

where 

f (K n ; k) = [t- 1 - G(k + K n ) - 5G(K n )]-\ (50) 
We also write the response functions Lfj. u in terms of cluster lattice Fourier transforms, i.e., 

L% v = -^[ dkLjjAV^-^. (51) 

On carrying out the sums over C and C in Eq. (|48p we thus obtain 

W°)=E [nfcS Knj K n , / rfkf(K n ;k)e i ( K " +lt )-^- R -)f(K n ,;k)e i ( K "' +lt )-^- R ^ 

—T T \JkL;u 

+ V /^I^Ekk, / t*f(K n ;k)e i < K » +i >^- R ^f(K n ,;k)e i ^' +£ )-( Ri - 11 ^ 

^ »r ™ "'fit 

-f^f iJ ]L A/A/;i ,(0) (52) 



K,L,M,N 



Since Lu ; „(0) is translationally invariant, L/j^(0) = Lj + / 1+ c :li j + / 1+ c' 1 ; l /(0) for a translation 
by an arbitrary lattice vector Rj = + R^ and we can sum it over the N c tile lattice 
vectors, L IJu (0) = jj-J2h Li+i u j+i i;v . Using this manipulation and applying the NLCPA 
coarse graining again so that e lkR/ ~ 1 we find 

W(0) = Ex rW ^+ E ^' L ' M 'V^' L ' J W;,(0) (53) 

K,L K,L,M,N 

with the x involving a convolution integral over the Brillouin zone, 



— E / dkf (K n ; k)e iK "- (R ^ Rx) f(K n ; k)e iK "-( R ^ Rj ) - f IK f LJ . (54) 
rz T , Jn t 



Ljj(0) is extracted by inverting a 'super' matrix, [1 — wx] which has dimension N c x N c x 
A/a x A/a, (A^a specifying the number of angular momentum quantum numbers). 

By comparing Eqs. ( 1461) and (1471) and using the definition L/j.,, = 0, the second term on 
the RHS of Eq. (1461) vanishes. We thus find the inter-cluster contribution to the conductivity 
to be 
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4m 2 / _ \UKL 

°l» = S S JlimMw t 1 - x(°H ^(O) JjKMMw ( 55 ) 

71-/2 " I,J,K,LAiMMM V /Ai,A 3 ,A4,A 2 

and the intra-cluster component to be 

= ^i,3 n ^2 P ~<c ^7,A4,Ai; M ^Ai,A2^KL,A2,A 3 ;^ 7 A3,A4- ( 56 ) 

7IYt \t 7c I,K,L Ai,A2,A3,A4 

VII. COMPARISON WITH THE CPA - ALLOYS WITH NO SHORT-RANGE 
ORDER. 

We have implemented the formalism outlined above using the Munich self-consistent, spin 
polarised, relativistic KKR (SPRKKR) code of Ebert et alM. Throughout we use an angular 
momentum cutoff l m ax = 3 which is necessary for studies of transition metal systems with 
significant d-electron weight in the electronic structure close to the Fermi energy. The current 
matrix elements which occur in the conductivity expression have odd parity and couple for 
example d-states to both p- and f-states. Omission of the effect of the latter can lead 
to an underestimate of the conductivity— Although in principle the recently developed 
SCF-KKR-NLCPA method^ can provide the appropriate self-consistent one-electron charge 
densities and potentials, p 7c (r/)'s and u 7c (rj) for our transport calculations, in these first 
applications we use those generated by the faster, simpler SCF-KKR-CPA method^ in 
order to explore the new aspects of our theory. 

Owing to the deviation of its resistivity from the common Nordheim (p ~ c(l — c)) 
behavior — and to earlier extensive studies made of itr^^^, the AgPd series of solid solutions 
provides the ideal initial test bed for our new method. We implement the formalism with 
NLCPA clusters containing 4 sites (Nc = 4) for this f.c.c. based system. Figure 1 shows 
the calculations of the resistivities of randomly disordered, substitutional Ag c Pd\^ c alloys 
compared with those calculated using the established CPA formulism of Butler—. Both 
implementations are presented with and without the so called 'vertex corrections'. We 
see that there is little difference between the two sets of calculations. Each describes the 
experimental trends well and the vertex corrections are found to be fairly insignificant for 
each approach. Both of these aspects have been discussed fully in relation to the CPA 
calculations in earlier publications (e.g.— so will not be discussed further here. Extensive 
work over 2 decades using the CPA formalism has shown how the resistivities of randomly 
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disordered transition metal alloys can be reliably described. The good agreement between 
our NLCPA results and those from the well established CPA method for these alloys where 
no short-range order is present is a very satisfactory first test of the new formalism. 




c 

FIG. 1: The resistivity of randomly disordered alloys Ag c Pdi- c alloys as a function of concen- 
tration c. The full lines are the NLCPA results (pink (green) lines with squares (diamonds) show 
those without (with) vertex corrections). The dashed curves show the CPA results (long dashes 
(dots) - without (with) vertex corrections). The experimental results of Guenault^ are shown for 
comparison (full blue lines with filled circles). 



VIII. THE EFFECTS OF SRO ON RESISTIVITY. 

Many properties of alloys such as resistivity are affected by short-range order. Indeed 
resistivity measurements are often used to monitor the changes in SRO which occur in an- 
nealing processes. If an alloy undergoes defect annealing after having been cold-worked there 
are significant changes in its physical properties owing to microstructural changes. For tech- 
nical applications it is important to know what these changes are so that physical properties 
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can be controlled. SRO plays an important role in this and resistivity measurements are 
used to follow its kinetics.— Our formalism is designed to help the interpretation of such 
measurements since it can describe the effects of short-range order on transport properties 
of alloys. It enables the calculation of the resistivity of a system to be made for a prescribed 
degree of SRO via the setting of the cluster configurational probabilities P 7C . Hence it can 
aid the extraction of SRO attributes from resistivity measurements. 

Our first application is to the b.c.c. based series of disordered Cu c Zn\- c alloys. We im- 
plement the NLCPA resistivity formalism using the smallest clusters and coarsest Brillouin 
zone tiling, i.e. Nc = 2. This means correlations only between nearest neighbors can be 
described. We incorporate SRO according to the following 3 prescriptions for the 4 config- 
urational weights, Pi = P(CuCu), Pi = P(CuZn), P3 = P(ZnCu) and P4 = P(ZnZn): 

• No SRO, P(CuCu) = c 2 , P(CuZn) = P(ZnCu) = c(l - c), P(ZnZn) = (1 - c) 2 . 

• Short-range order (minimising number of like nearest neighbors), for c > 0.5, 
P(CuCu) = (2c - 1), P(CuZn) = P(ZnCu) = (1 - c), P(ZnZn) = and for 
c < 0.5, P{CuCu) = 0, P(CuZn) = P(ZnCu) = c, P(ZnZn) = (1 - 2c). 

• Short-range clustering (maximising number of like nearest neighbors), P{CuCu) = c, 
P(CuZn) = P(ZnCu) = 0, P(ZnZn) = (1 - c). 

Figure 2 summarises our findings. Vertex corrections are included and are shown to be large 
for these systems. This concurs with earlier results for the randomly disordered alloys which 
have low d-electron weight in electronic structure around the Fermi energy. In the absence of 
any short-range ordering, the results show approximate adherence to the expected Nordheim 
c(l — c) behavior and the results both with and without vertex corrections are very close 
to the CPA results and experimental results 60 . Incorporating short-range order with extent 
only between nearest neighbors decreases the resistivity as expected for all concentrations. 
The resistivity now follows a rough c(l — c) + Ac 2 (l — c) 2 dependence so that the greatest 
reduction occurs at the stoichiometric concentration of c = 0.5. Conversely when short- 
ranged clustering is included the resistivity increases for all concentrations and the c(l — c) 
behavior returns. Evidently current is enhanced on both types of atom when they are 
surrounded by unlike neighbours. These results are in good agreement with experimental 
measurements of the resistivity which show the resistivity to decrease significantly when the 
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alloys are annealed so that short-ranged order is induced. 60 




FIG. 2: The resistivity of Cu c Zn\- c alloys as a function of c, concentration. The full (red) line 
with crosses shows the calculations when no SRO is included. (These are nearly indistinguishable 
from the CPA results (green line)). The (green) line with asterisks shows the NLCPA results when 
SRO is included and the square boxes (light blue line) show results when short-ranged clustering is 
included. Results are also shown where the vertex corrections have not been included: full (yellow) 
line with plus signs - NLCPA results for no SRO, (the dashed line shows the CPA results); lines 
with diamonds - NLCPA results including SRO; and lines with triangles - NLCPA results including 
short-ranged clustering. 

We have also investigated the effect of short-range order on Ag c Pd\- c alloys by once more 
choosing configurational weights such that the number of like neighbours is minimised for 
each concentration. Figure 3 contains the results. The effect of SRO is less than that found 
in CuZn. Below c = 0.3 the effect is negligible whereas for larger concentrations short-range 
order depresses the resistivity a little showing how the current is enhanced on a site when 
surrounded by unlike neighbors. Experimental measurements on Ag c Pd\_ c alloys find that 
annealing has a smaller effect 587 than in Cu c Zn\- c in line with our calculations. It is also 
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found that cold work causes little change to the resistivity suggesting that additional defects 
such as dislocations may already be present affecting the measurements. This may indicate 
the origin of our underestimate of the resistivity shown in Figure 1 and Figure 3. 




c 

FIG. 3: The resistivity of Ag c Pd\- c alloys as a function of c, concentration. The full (green) 
line with diamonds shows the calculations when no SRO is included. The (red) line with triangles 
shows those when SRO is included. Vertex corrections are included in both plots. The experimental 
data 58 are also shown (blue line with filled circles). 

IX. CONCLUSIONS 

Short-range ordering and clustering dramatically affect the transport properties of many 
alloys. Indeed resistivity measurements, which can be made easily and rapidly, provide 
a good way to monitor microstructural changes that occur in materials processing. In 
this paper we have described a way to make quantitative calculations of the resistivity 
of disordered systems which possess short-range order or clustering. The ab-initio theory 
starts from the density functional theory for these systems devised recently by Rowlands 
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et al.— using the SCF-KKR-NLCPA electronic structure method. Our first calculations for 
Cu c Zrii- c and Ag c Pd\_ c show the expected decrease of resistivity when short-range order is 
imposed whereas short-ranged clustering produces an increase. For the randomly disordered 
alloys the results are very similar to those that have been produced from established KKR- 
CPA resistivity calculations based on Butler et al. work— >2i. These first calculations have 
not included the effects of the short range order on the self-consistent charge densities and 
potentials that are available from the SCF-KKR-NLCPA method. So far also short-range 
clustering and ordering effects over only the shortest nearest atomic neighbor range have 
been included. The computational development work is in progress to remove these current 
practical limitations. The counterintuitive behavior of the resistivities of the K-state alloys 
such as NiCr,NiMo and PdW will be ideal next systems to study. 
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